Conformational transformations induced by the charge-curvature interaction 
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A simple phenomenological model for describing the conformational dynamics of biological macro- 
molecules via the nonlinearity-induced instabilities is proposed. It is shown that the interaction 
between charges and bending degrees of freedom of closed molecular aggregates may act as drivers 
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giving impetus to conformational dynamics of biopolymers. It is demonstrated that initially circular 
aggregates may undergo transformation to polygonal shapes and possible application to aggregates 
of bacteriochlorophyl a molecules is considered. 

PACS numbers: 



I. INTRODUCTION 



The functioning of biological macromolecules is determined by their tertiary structure and different types of semi- 
flexible polymers in solutions exhibit a variety of conformational phase transitions 1\. Even modest conformational 
changes modify long-range electronic interactions in oligopeptides j^, they may remove steric hindrances and open 
' the pathways for molecular motions which are not available in rigid proteins 3]. In particular, it has been recently 

o ' 



shown m that flexibility increases the hydrogen accessibility of DNA fragments and in this way facilitates strand 
breaks in DNA molecules. There is also a strong belief that the conductivity of DNA is due to thermal motion of 
£3 ' small polarons 0, IE Q • 

i-^j ■ Research in solitonic properties of the chains with a bending has been initiated in the recent years 0, 0, lldl 111! Il2| . 

In particular, it was shown that the bending of the chain could manifest itself as an effective trap for nonlinear 
§ i excitations gQQ and that the energy of excitations decreases when the curvature of the bending increases Q 
^ , Quite recently, there has been a growing interest in studying nonlinear charge and energy transport in soft condensed 

systems (polymers, membranes) with self-consistent account of coupling between nonlinear excitations and the shape of 



the systems. It has been demonstrated |l8j that a mismatch of length scales in the presence of magnetic solitons leads 
to an elastic deformation on a soft magnetic surface. A phenomenological model for describing the conformational 
dynamics of biopolymers via the nonlinearity-induced buckling and collapse instability was proposed in |2fJ | . It was 
shown there that the nonlinear excitations may cause local softening of polymer bonds. That is the effective bending 
rigidity of a chain may become negative nearby the nonlinear excitation and in this way the nonlinear excitation 
causes a buckling instability of the chain. 

In this paper we study a simple model for electron-curvature interactions on closed molecular aggregates. In 
particular we show that due to the interaction between electrons and the bending degrees of freedom the circular 
shape of the aggregate may be become unstable and the aggregate takes the shape of an ellipse or a polygon. It is 
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shown that the interaction between complexes may stabilize the ringlike shape of the The paper is organized as follows. 
The theory is applied to study the conformational transitions in light-harvesting complexes in purple bacteria. In Sec. 
II we describe a model. In Sec. Ill we present an analytical solution to the Euler-Lagrange equation. In Sec. Ill we 
compare our analytical results to results obtained directly by numerical simulations. In Sec. IV the interplay between 
shape of complexes and intercomplex interaction is studied. Sec. V is devoted to application to conformational 
transformations in light-harvesting complexes. Sec. VI presents some concluding remarks. 



II. THE MODEL 

Let us consider a chain consisting of L units labelled by an index n, and located at the points r n — {x n , y n , z n }. 
We are interested in the case when the chain is closed and so we impose the periodicity condition on the coordinates 

r„ 

f n = f n+L . (1) 

The chain flexibility is accounted for by employing a microstructure consisting of many sequentially joined rigid rods 
and by incorporating a bend potential at each point of rotation [2 lj | . Thus the Hamiltonian of such a polymer chain 
has the form 

H = U + H eh (2) 
with the potential energy of inter-unit interactions U(...,r n , r n +i, ■■) which we take in the form 

U = U S + U B (3) 

where 

Us = \ ]T(|r„-r; +1 |-a) 2 ( 4 ) 

n 

is the stretching energy in the harmonic approximation. Here a is an equilibrium distance between units (in what 
follows we assume a = 1) and a is a dimensionless elastic modulus of the stretching rigidity of the chain, and 

^ = ^ 1- (5) 

n "'n/ ""max 

is the bending energy. Furthermore, 

K n = I •?>, - r„_i I = 2 sin ~y (6) 

determines the curvature of the chain at the point n. The vector 

- _ r n+ i - r n 

TVl I -* _ I 

\r n+ i - r n \ 



(7) 



is the unit tangent vector, a n is the angle between the tangent vectors f n and r„_i, n m ax — 2 sin (a max /2) with a max 
being the maximum bending angle, k is the elastic modulus of the bending rigidity of the chain. 

We assume that there is a small amount of extra electrons ( or holes) on the chain. The Hamiltonian of electrons 
and their interaction with conformational degrees of freedom is 

H e i = Hhop + H e i- con f, (8) 
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where 



Hhop J 



\lp n - 1p n+ l 



describes the motion of electrons along the chain, and 

1 

2 ^ 



H, 



el—conf 



(IV'n+ll 2 



(9) 



(10) 



gives the interaction between electrons and bending degrees of freedom. In Eqs. (|8|)- (|1U|) rp n {t) is the electron wave 
function at the site n, the parameter J describes the electron hopping in the chain, and x is the curvature-electron 
coupling constant (see Appendix for details). 
The quantity 



(11) 



gives the total density of extra electrons which can move along the chain and participate in the formation of the 
conformational state of the system. We will neglect the interaction between electrons. It is legitimate when the total 
density of electrons in the chain v is small. Combining Eqs. JSJ and (|10fl . we notice that the effective bending rigidity 
changes close the points where the electron is localised. For positive values of the coupling constant \ there is a local 
softening of the chain, while for x negative there is a local hardening of the chain. 

In what follows we assume that the chain is planar (z n = 0) and inextensible (<r — > oo): 

\r n -r n+1 \ = l. (12) 
III. CONTINUUM APPROACH 



We are interested here in the case when the characteristic size of the excitation is much larger than the lattice 
spacing. This permits us to replace ip n (t) by the function ip(s,t) of the arclength s which is the continuum analogue 
of n. Using the Euler-Mclaurin summation formula we get 



H = U b + H e 



(13) 



Ub = — k J k(s) 2 ds 



(14) 



H p1 = < J 



d s ip 



X K{s) 2 \iP\ 2 \ds 



(15) 



Being interested here in small curvature effects: k(s) <^ n max , we consider the bending energy in the harmonic 
approximation given by Eq. (|14f> . 



A. Ground state of the chain 

The electron ground state wave function (£>(s) and the shape of the chain f(s) may be obtained by minimizing the 
functional 



(16) 
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v= - I (/){s) 2 ds, (17) 



with H e i and Ub given by Eqs l|15|) and 1(14(1 under the constraint 

L 

L 

o 

which is a continuum analog of Eq. . The inextensibility constraint 1)12(1 reads in the continuum limit 

\d s rf = 1. (18) 
The inextensibility constraint 118(1 is automatically taken into account by choosing the parametrization 

d s x{s) — sin 9(s), d s y(s) — cos 8(s) (19) 
where the angle 9(s) satisfies the conditions 

6{s + L) = 2tt + 6{s) (20) 

and 

L L 

J cos 6{s) ds = J sin 8{s) ds = (21) 



which follow from Eq. Note that, in the continuum limit, the curvature of the chain k(s) given by Eq. © can 
be expressed as 

k(s) = |d s V( S )|. (22) 
Thus, in the frame of the parametrization l(19|) . k(s) = d s 9 and the functional ((16(1 takes the form 



£ = J\jv (ds^f + Q - X *V) (9 S ^) 2 U, (23) 
where the rescaled function </?(s) = \pv (j)(s) which satisfies the normalization condition 

L 

\ J ^(*) = 1, (24) 
o 

has been introduced. 

The Euler-Lagrange equations for the problem of minimizing £, given by Eq. ((23(1 under the constraint ((24(1 become 

d 2 sV > + j {d s 0f <j) - Xip = 0, (25) 

d s (d s e (l-wtp 2 )) =0 (26) 

where A is the Lagrange multiplier and 

w = —— (27) 
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is a coupling constant which characterizes the strength of the charge-curvature interaction in terms of the bending 
rigidity of the chain and the charge density. We are interested in solutions of Eq. (|25(l subject to the periodic boundary 
conditions 

p(a) = + (28) 
n 

where n is an integer which characterizes the shape of the chain (see below) . Integrating Eq. H26[l , we get 

A_ 

1 — w tp 2 

where A is an integration constant. Taking into account the condition l|2l)fl we obtain that the integration constant 
A is determined by the relation 

A=^I (30) 



d s0 = ; —2 ( 29 ) 



where the functional / is given by the relation 



L 

1 1 f ds 



I L J 1 — w <p 2 
o 

From Eqs. (|19() and 129(1 we see that the shape of the chain is determined by the equations 

s s 




(31) 



^ s ) = 77 / 1 1 -TT^ ds '- ( 32 ) 

LI J 1 — w {p z {s') 







B. Solution of the Euler-Lagrange equations 

There are two kinds of solutions to Eqs (|25|l and (|29|) . 

• Circular chain. 
Charge is uniformly distributed along the chain 

tp = 1 (33) 
where the normalization condition l|17l) has been used, and the curvature of the chain is constant 

k(s) =d s 6= —— . (34) 
1 — w 

This case corresponds to a circular chain 

s s 
x = R sin — , y = —it cos — . 

R' y R 

The radius R of the circle can be obtained by puting Eq. into the boundary condition l(2l7|) . As a result we 
have 

R=±-. (35) 

Z7T 



The energy of the circular chain is thus 



2tt 2 

£circ — —j— k(l — w). (36) 
Li 
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• Polygonally deformed chain. 
Let us consider now the case of spatially non- uniform distributed electrons. Inserting Eqs. 1)29(1 - 1131(1 into Eq. 
(EU we get 



£ = J v J (d s ip) 2 ds + —j-— I- 



(37) 



We restrict our analytical consideration to the case when the charge-curvature coupling is weak and/or the 
charge density is low: w <C 1. Expanding the functional / in terms of the small parameter w we obtain from 
Eq. 



£ = 



2ir 2 k 
L(l + w) 



where 



G 



2X 2 



(38) 



(39) 



Jk(l +w) 2 

is an effective nonlinear parameter. For small w one can neglect the last term in Eq. I(38|) and the Euler-Lagrange 
equation for the functional ((38(1 then takes the form 

Gv 



9 s V + 2— ^ 3 -A^-0. 



Straightforward calculations show that Eq. ((40(1 has a solution of the form 

/ ^ J 
<p = R\l — ; —. an 



A 



(40) 



(41) 



{2-m)Gv \ V (2-m) 

where dn(u\m) is the Jacobi elliptic function with the modulus m [2^. Inserting Eq. ((41|) into the boundary 
condition 1(28(1 and the normalization condition 1(24(1 . we find that the Lagrange multiplier A and the modulus 
m are determined by the equations 



A L 



(2 — m) n 



2K(m) 



(42) 



Gv = —K(m)E(m) 



and the charge distribution along the chain is given by 



K (2nK 
E dn \— S 



(43) 



(44) 



where K(m) and E(m) are the complete elliptic integrals of the first kind and the second kind, respectively 2 
The curvature of the chain is given by the equation 



, s r, „ 1 ( K T 2 (2nsK 
k{s) = d s 9^ — I [l + w — dn 2 
R \ E 



L 



Integrating Eq. J25J, we get 



e( S )=^l(s + W J^E(a\m) 



(45) 



(46) 



7 



where E (aim) is the incomplete elliptic integral of the second kind, and a — am ( 2n L K s|m) is the amplitude 
function |22j ■ By using the Fourier expansion for the amplitude function for small m we obtain from Eq. 146(1 



fi(s) s; — (l+ui)/s + — /msinl — si . (47) 

Inserting Eq. 147(1 into the closure condition 121(1 . we find that it is satisfied for n > 2. Eqs. ((32(1 and 147(1 
describe a polygon: for n = 2 it is an clliptically deformed chain, while for n = 3 it has a triangular shape (see 
Fig. HI). We see from Eqs. ^4*4*jl and lf4"5|) that the polygon structure is a result of the self-consistent interaction 
between electrons and bending degrees of freedom: extrema of the curvature and of the charge density correlate: 
in the case of the softening electron-curvature interaction (x > maxima of curvature and charge density 
coincide, while in the case of the hardening interaction the minima of the curvature coincide with the maxima 
of the charge density. Eq. ((43(1 shows that, for a given value of the nonlinear parameter G, the n-gon structure 
appears when the charge density exceeds the threshold value v n : 

n 2 

V > Vn = 4G (48) 



The energy difference between the n-gon structure and the circular chain is given by the expression 

4tt 2 GJi 
~3L ~E i 



4_2 Q j 2 

£n - £ cl rc = -z-r (3E 2 - (2 - m)EK - (1 - m)K 2 ) . (49) 



The normalized energy difference 



F — F 

An = — -p (50) 

&circ 



for n — 2,3 versus the charge density is shown in Fig. We note that when the charge density is above 
the critical value the deformed structure with spatially inhomogeneous charge distribution is energetically more 
favorable than the circular system with a uniformly distributed charge. The state with elliptically deformed 
chain n = 2 is the energetically most preferable. 

Note also that our analytical approach was based on the assumption that wip 2 <C 1. Taking into account Eq. 
(gU), this means that it is legitimate to consider not too sharp distributions which correspond to wK(m) <C 1 
or m <C 1 — exp{— 0.72/iu}. 

IV. NUMERICAL STUDIES 

To check our results we have performed also several numerical studies. To this end we carried out the dynamical 
simulations of the equations 

d _ dH 
dt or n 

l Jt^ = -dr n (52) 

with the Hamiltonian H being defined by Eqs. l(2|l- l(l(J|) . Thus the conformational dynamics is considered in an 
overdamped regime with the friction coefficient r\. Then we took as our starting configurations systems involving the 
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electric charge density of (almost) the same magnitude {ip n ) at all points (we broke the symmetry by increasing the 
density at one point of the chain by 1%). Initially, all the lattice points were placed at symmmetric points on the 
circle of an appropriate radius (see Fig|3J). We performed such simulations for several values of the charge density. 
Due to the absorption the energy of the system was decreasing during the evolution and the system was evolving 
towards a minimum. At the same the points of the chain were moving from their initial to their new positions. 

We considered both the cases of the hardening and of the softening electron-curvature interaction. 

In the case of a hardening electron-curvature interaction (x < 0) a typical final distribution of the chain points 
is shown in Fig^l The charge density and the curvature distributions are in full agreement with the results of our 
analytical considerations: the curve is more flat where the density of the electrons is maximal. In fact, the ellipse-like 
shape is rather robust as it arises for a large range of parameters (of the strength of the hardening electron-curvature 
interaction and of the anharmonicity coefficient n max ). 

Complexes with a softening electron-curvature interaction are much more flexible. Their equilibrium shape depends 
drastically both on the anharmonicity and on the charge density. Figs [5] and [B] demonstrate how drastically the shape 
of the complex and the charge distribution along the chain can change as a function of the total charge density v: 
increasing the total charge by 5% can lead to the localisation of almost the whole charge of the system at one place. 

In our numerical work we also studied the stability of our 'final' field configurations - i.e. the configurations which 
we thought the system was settling at. This we studied by perturbing the system. Such perturbations were introduced 
in two stages. First we changed the electric charge of the configuration by multiplying all 'final' values of the electric 
charge by a constant factor /i; this had the effect of changing the energy of the system. Then we performed the new 
minimisation and, when the system appeared to have settled at the new 'final' configuration, we changed back its 
ip n by a new multiplication by As |VVi| 2 is conserved during the evolution, the final system had the same 

value of it as the original 'unperturbed' fields. The results of the further minimisation were then compared with the 
original 'final' fields. 

When we applied this technique to our field configuration shown in Fig. we found that the system was really 
unchanged by this perturbation; in fact the perturbation led to an overall rotation of the system by one lattice point, 
but the sequence of values of the fields was essentially the same thus showing the stability of the found minimum. 



The aim of this section is to investigate how the interaction between complexes influences the shape. We will 
consider the system which is described by the Hamiltonian 



where £j is the energy of the j-th aggregate which is given by Eq. (12MI) with tp replaced by fj and 9 replaced by 0j, 
and Uij is the interaction energy between particles. The latter we will take in the form of the Gay-Berne potential 
[2^ which is a generalization of the Lennard- Jones 12-6 potential and is widely used to study translational and 
orientational ordering in systems of aspherical molecules. We consider small deviations from the ringlike structure of 
aggregates and so we neglect the difference of the well depths for side-to-side and end-to-end configurations. In this 
case the Gay-Berne potential between two parallel uniaxial molecules is given by 



V. EFFECTS OF INTERCOMPLEX INTERACTION 




(53) 
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where ru = fijf^j is the interparticle vector, 



<t [fijj = ctq [1 - j— [fij ■ ej (55) 



i + c 

is the anisotropy parameter where e is a unit vector specifying the axes of symmetry. The anisotropy coefficient ( is 
determined by the lengths of the major and minor axes <7|| and a± 

T 2 



C = %~^2 (56) 



a \\ 

and Co gives a charasteristic length scale while Uq determines the intensity of the interaction. 

The centers of densely packed circular aggregates of the radius R create a two-dimensional triangular lattice fj = 
ji^i +i2«2, (j = {jx-,32), 31,3-2 = 0, ±1,±2, ...) with the basic vectors a\ = £ (1,0) and a 2 — £ (1/2, V3/2) where 
£ is the lattice constant (see Fig. EJ- Being elliptically deformed in such a way that the major axes of all aggregates 
are parallel to the cc-axis, the centres of densely packed elliptical aggregates create a lattice fj = j\b\ with the 

basic vectors b\ = £ (1 + u, 0) and b 2 = I (1/2 (1 + u), s/3/2 (1 — it)) (see Fig. [SJ, where the parameter u is given by 

U =^-^. (57) 

Now we study the ground state of this system by using a trial function approach. We assume that this state is 
spatially homogeneous and relying on the results of the previous section, we assume that the electron trial function 
and the trial curvature can be taken in the form 

ifj = ^cosa + \pi sin a cos (— j j , (58) 



^, = I(l + 7 cos(|)) (59) 

The functions (|58[) and (j59(l can be considered as a truncated Fourier expansion of the solutions j|41|) and 145|l in which 
the coefficients a and 7 are variational parameters and R is the radius of the cylindrically symmetric aggregate. The 
function (|58|l satisfies both the periodicity condition (|28|l and the number of particles constraint Ijl7|l . In the limit 
7 < 1 the shape of the curve with the curvature given by Eq. I|59|l is parametrically determined by the expressions 



*>-*((! + J> <»(£)+ £«»(£)) 



„(.)=*((!- 2) + (if)). (60) 

Thus the lengths of the major and minor axes of the curve H60JI are given by 

a ]} =R(1 + 1), a ± = R(l-l) (61) 

and comparing Eqs (|57|) and l|t)l|) . we see that u = 7/3. 

Inserting Eqs. 1)58(1 and (|59|l into Eqs. ((53(1 . 1(23(1 and 1(54(1 for an energy per aggregate we get 

w = Sir + U (62) 

1* a 

where 

£ tr = ^-l 1 + -J- - ±£v (8 + 5 7 2 - 7 2 cos(2a) + 8V2 sin(2a)) + ^- v sin 2 a \ (63) 
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is the energy of an isolated aggregate, and 

5 

U = Y / U(A 3 ) (64) 

is the energy due to the interaction between aggregates in the lattice. In Eq. i|64|) the function U (Aj) is given by Eqs. 
1)540. (|55|l with e = (1, 0) and vectors Aj = <tq ((1 + 3) cos (^J , (1 — ^) sin (^j) connect nearest and next-nearest 
neighbours of the lattice. The interaction energy l|64(l has a minimum at i — 2 1 / 6 7 = which corresponds to a 
system of densely packed circular aggregates. Expanding the function (|64[1 in the vicinity of this point, in powers of 
the variational parameter 7, we get 

U = -^U + cU j 2 + ■■■ (65) 
where the numerical coefficient c 2.05. According to the variational principle we should satisfy the equations 

d a £tr = 0, (66) 



a 7 {£ tr + U)=0. (67) 

From Eq. (|B1)|) we get 

Inserting Eq. I|68|l into Eq. I|63|) and expanding it in terms of 7 we obtain that 

£tr=£ cl rc + ^(i-tt) 7 2 +-B 7 4 + --- (69) 

where 



2R V v c 



2x J + X 



' (70) 



is the critical charge density and the notation 



16 i? J 3 ' 

is introduced. Thus in the framework of the variational approach the energy of a single complex has a single minimum 
at 7 = when v < v cr and in this case the aggregate has a ring-like shape. When v > v cr and 2 x > J the energy 
(16911 possesses two equivalent minima with ±70 (70 7^ 0) ( see Fig. |§J. As it is seen from Eqs. (|60|l the finite value of 
70 implies that the aggregate is elliptically deformed either along the a;-axis (when 7 > 0) or along the y-axis (when 
7 < 0). Note that in the limit x ^* J> v cr coincides with j/ 2 given by Eq. (|4*51) . Combining l|6"5|l and l|6^l) we see that 
the interaction between aggregates modifies the condition for appearance of the low-symmetry form. Indeed, even 
in the case when v > v cr (an isolated aggregate has an ellipse-like shape) in the condensed phase of aggregates for 
strong enough inter-aggregate interaction we may have an inequality 

u< ( 1 + ^JT R ) (71) 

which means that interacting aggregates are ring-like and create a densely packed crystallic structure with the group 
symmetry D 6 h (see Fig. [7||.When 
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the ellipse-like shape survives in the condensed phase which has a rectangular unit cell and transforms in accordance 
with the group symmetry Dih ( see FiglHJ . There are two equivalent types of arrangement in the low-symmetry phase 
when the aggregates are elliptically deformed either along the £-axis (when 7 > 0) or along the y-axis (when 7 < 0). 
In accordance with this the condensed phase of the aggregates must have a domain structure, i.e. it must consist of 
various regions in which the direction of long axes are different. 



VI. APPLICATION TO BACTERIOCHLOROPHYLL a COMPLEXES 



The X-ray crystallography shows 23] that bacteriochlorophyll a (Bchl a) molecules in the LH complex are organized 
in two concentric rings: the B800 and B850 rings. The former consists of nine well-separated Bchl a molecules with an 
absorption band at 800 nm and the latter consists of 18 Bchl molecules with an absorption band at 860 nm (see Fig. 
1 in Ref. Measurements of the anisotropic properties of the absorption of isolated LH2 complexes from bound to 

mica surfaces Il3l . the fluorescence-excitation spectra of individual LH2 complexes from bacteria Rhodopseudomonas 
acidophila 

QE30 

showed that the complexes are generally not cylindrically symmetric but reveal a deformation 
of the circular complex into a shape with C2 symmetry. Quite recently the shape of the LH2 complex from bacteria 
Rhodobacter spheroides 2.4-1 in detergent solution has been determined from synchrotron small-angle X-ray scattering 
data [l?! ]. It was shown that, in contrast to the cylindrical crystal structure with a diameter of 6.8nm, the shape of 
an isolated LH2 complex is an oblate plate with an eccentricity e = 0.59. It was conjectured in Ref. Q| that the 
extremely dense packing of LH2 in the crystals causes the cylindrical symmetry of the complexes. A variety in shapes 
and conformations for light-harvesting LH1 complexes and different types of their packing in two-dimensional crystals 
revealed by atomic force microscopy were reported recently in Ref. |24| . 

A qualitative explanation of this phenomena follows from the theory developed in previous sections. It has shown 
that in the presence of sufficiently strong charge-curvature interaction an isolated complex has an ellipse-like shape 
while the interaction between complexes in the form of an anisotropic Gay-Berne potential stabilizes the ringlike 
shapes of the complexes. However if the intensity of the intercomplex interaction is less than some threshold value 
(see Eqs. and (Hp) the non-circular shape of the complex is preserved in the condensed phase. 

The model is too simplified for quantitative predictions for light-harvesting complexes but we believe that it contains 
interesting physics which should be important for further studies of such systems. 



VII. CONCLUSIONS 



In this paper, we have investigated the role of the electron-curvature interaction on the formation of the ground 
state of closed semiflexible molecular chains. We have found that the coupling between electrons and the bending 
degrees of freedom of the chain can induce a local softening or hardening of chain bonds, i.e. the effective bending 
rigidity of the semiflexible chain changes as the density of charge changes along the chain. When the charge density 
and/or the strength of the electron-curvature coupling exceed a threshold value, the spatially uniform distribution of 
the charge along the chain and the circular, cylindrically symmetric shape of the chain become unstable. In this case 
the ground state of the system is characterized by a spatially non-uniform distribution of electrons along the chain 
and the chain takes on an ellipse-like (or in general polygon-like) form. 
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APPENDIX A 

The electron-conformation interaction can be obtained from the Coulomb interaction between an electron and 
charged groups of the molecule 

H int = J2 v (\r n -f n/ \) |V™ I 2 (Al) 

n,n f 

where the matrix element V (|r„ — r n > |) describes the Coulomb interaction between an electron which occupies the 
site n in the chain and the charged group at the site n' . From the inextensibility condition (|12f) we see that the nearest 
neighbour interaction terms (|n — n'\ = 1) do not contribute to the electron-conformational interaction. Taking into 
account that 

{r n +2 - r n ) = 2a 2 + 2 (f n+1 - r n ) ■ (f n+2 - r n +\) = 4a 2 - a 2 n 2 l+1 (A2) 

we obtain from Eq. (|A1|) that, in the next-nearest-neighbour approximation, the energy of the electron-conformational 
interaction takes the form 

Hel-conf-app = ~ ~ ^ X + l^nP (A3) 

n 

where the notation 



1 dV 

X = 

has been used. 



(A4) 

r=2a 
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FIG. 2: The normalized energy difference A„ from Eq. l|5()l for the ellipse-like n = 2 ( solid line ) and triagon-like chain n — 3 
(dashed line) with 2\ — 4 J = k. 
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FIG. 3: The top panel shows the initial shape of the chain; the bottom panel shows the initial charge (solid line) and curvature 
(dashed line) distribution along the chain. 
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FIG. 4: The top panel shows the equilibrium shape of the chain, and the bottom panel shows the charge distribution 
(solid line) and curvature variation (dashed line) along the chain in the case of hardening electron-curvature interaction with 
v = 0.22, X = -2, k = 1J = 0.4 
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FIG. 5: The top panel shows the equilibrium shape of the chain and the bottom panel shows the charge distribution (solid 
line) and curvature variation (dashed line) along the chain in the case of softening electron-curvature interaction with v — 
0.22, x = 2.15, k = J=l 




FIG. 6: Same as Fig. CD with v = 0.23, % = 2.15, k = J = 
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FIG. 7: Arrangement of densely packed ring-like aggregates. 




FIG. 8: Arrangement of densely packed ellipse-like aggregates. 
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FIG. 9: Energy per aggregate (Eq. 16210 as a function of the trial anisotropy parameter 7 for the charge density below the 
threshold (dashed) and above the threshold (solid). 




